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Abstract 

We have performed a search for radiative 6-hadron decays using events pro- 
duced in pp colUsions at ^/s = 1.8 TeV and collected by the Collider De- 
tector at Fermilab. The decays we considered were — K*^{-^ K~'k~^)^, 
— *■ K'^K~)'y, A° p7r~)7, and their charge conjugates. Two 

independent methods to identify photons from such decays were employed. 
In the first method, the photon was detected in the electromagnetic calorime- 
ter. In the second method, the photon was identified by an electron-positron 
pair produced through the external photon conversion before the tracking 
detector volume. By combining the two methods we obtain upper limits on 
the branching fractions for the D% and Aj^ radiative decays, which, at 
the 95% confidence level, are found to be 8(3^ K j) < 1.4 x 10"*, 

^ (^7) < 1.6 X 10-'', and H(Ag ^ A7) < 1.9 x 10"=^. 
PACS numbers: 14.40.Nd, 14.20.Mr 
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I. INTRODUCTION 



Flavor-changing neutral currents (FCNC's) are suppressed in the Standard Model (SM) 
by the Glashow-Iliopoulos-Maiani mechanism [Q, and such transitions can only result from 
higher order processes. The "penguin" process is one such example, where an effective FCNC 
b ^ s OT b —>■ d transition proceeds through the emission and reabsorption of a virtual W 
boson. A photon, gluon, or Z boson is emitted from the quark or the W in the loop, with 
the presence of a photon signaling an "electromagnetic" penguin process (see Figure |I|). 

It is expected in the SM that the top quark dominates in the fermion part of the loop 
of the diagram. The existence of non-SM heavy charged particles, however, could affect 
the branching fraction for this decay. In addition, direct CP-violating effects could be 
enhanced by processes beyond the Standard Model. Therefore, measurements of radiative b 
hadron decays, constitute low energy probes for physics beyond the SM p|. Within the SM 
framework, radiative b ^ s decays are sensitive to the magnitude of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element |Vts|, while radiative b ^ d decays are sensitive to \ Vtd\- 
Ratios of branching fractions involving b d'y and b ^ s'y decays can thus be used to 
measure the ratio (|Vjd/|Vts|)- This ratio determines the length of one side of the unitarity 
triangle, and may explain the source of CP violation in the SM 

The branching fraction for the exclusive radiative decay K*^'-/ was first measured 

by CLEO to be (4.55l[Ji^±0.34) x 10^^ [|. The most precise measurements of the branching 
fraction are (4.23 ± 0.40 ± 0.22) x 10~^ by the BABAR collaboration § and 

(4.96 ± 0.67 ± 0.45) x 10"^ by the BELLE collaboration 0. Both collaborations have also 
measured the branching fraction B{B^ K*"'-/); with (3.83 ± 0.62 ± 0.22) x 10"^ obtained 
by BABAR ||, and (3.89 ± 0.93 ± 0.41) x 10'^ obtained by BELLE §. BELLE has also 
reported B{B pl)/B{B K*-f) < 0.19 at 90% confidence level (CL) 0. The branching 
fraction for the inclusive radiative decays B X^'y, where Xg represents a collection of 
hadrons containing strange quarks, was also measured by CLEO to be (3.15 ± 0.35 ± 0.32 ± 
0.26) X 10~^ 0, where the first uncertainty is statistical, the second is systematic, and the 
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third is for model dependence. The studies of the heavier 6-hadron decays such as and Ab, 
which are not produced at the T{AS), must be done at the higher energy machines, such as 
the Tevatron. No exclusive radiative decays of nor A° have been observed to date. From 
a search for 5^ 07 decays, the DELPHI collaboration obtained 07) < 7.0 x 10~^ 

at 90% CL [|. 

Even though calculations for the exclusive decay rates have higher theoretical uncertain- 
ties compared to inclusive decay rates, ratios of exclusive b ^ d'j and b ^ S'y branching 
fractions can be calculated with good precision and the determination of (| Vf(i/|\4s|) is fea- 
sible with the use of exclusive decays [0. This is especially useful for a hadron collider 
environment, where the experimental signature for radiative b decays is much cleaner when 
exclusive decays are considered. 

In this paper we report the results of a search for ^ K {-^ K~7r'^)'j, B^ <t>{~^ 
K^K~)'y, and A° — > A(— > pn^)'y decays in events produced in pp collisions at y/s = 1.8 TeV 
and recorded by the Collider Detector at Fermilab (CDF) during 1994-96. Two methods 
to identify such decays are employed. In the first method (Method I) [|1^], the photon is 
detected in the electromagnetic calorimeter. The trigger for this method required a minimum 
energy deposition in the calorimeter and two oppositely charged tracks that were distinct 
from the calorimeter signal. In the second method (Method II) ||I2[, the photon is identified 
by an electron-positron pair produced through an external photon conversion within the 
tracking detector volume. One of the conversion electrons, detected in the electromagnetic 
calorimeter, served as a trigger for recording these events. The b hadrons are then exclusively 
reconstructed with four charged tracks. 



II. COLLIDER DETECTOR AT FERMILAB (CDF) 



Since CDF is described in detail elsewhere |jT3[, we describe here only the components 



relevant to this work. In this paper we use a cylindrical coordinate system (r, 0, z) with the 
origin at the nominal interaction point, the z axis parallel to the nominal beam direction, 
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r the distance from the beam in the plane transverse to the z axis, and the azimuthal 
angle. We define 9 to be the angle with respect to the +z direction and the pseudorapidity 
as r/ = -ln[tan(6'/2)]. 

The tracking systems consist of a silicon vertex detector (SVX), a vertex time projection 
chamber (VTX), and an open-cell multiwire drift chamber (CTC), all immersed in a 1.4 T 
solenoidal magnetic field aligned with the z axis. The SVX [l^ is the innermost system, 



with its four layers of single-sided silicon microstrip detectors in the radial range of 3.0 to 
7.9 cm. The active area is 51 cm long in z and covers 60% of the pp interaction region. 
The microstrips all run parallel to the z direction and therefore track charged particles in 
the transverse plane. The SVX measures the impact parameter of tracks with respect to 
the beam line with a resolution of (Jd{PT) = (13 + 40/pt) /wni, where pt is the momentum 
of the track in the transverse plane in GeV/c. This precision close to the beamline helps 
distinguish the tracks of B decay products from those originating at the pp interaction point. 



The VTX [0 surrounds the SVX and consists of 28 drift modules with an outer radius 
of 22 cm and z coverage up to ±136 cm. The VTX tracks particles in the r — z plane 
and provides a measurement of the actual pp interaction point along the z axis with a 
resolution of 1 to 2 mm. From a combination of this information with SVX measurements, 
the transverse beam profile has been measured with an accuracy of 25 fim. 

Outside the VTX lies the CTC |l^, which extends out to a radius of 138 cm and 
\z\ < 160 cm. It contains 6156 wires arranged in 84 layers, which are further grouped into 
9 "superlayers" . Five of these superlayers are made of twelve layers of wires strung parallel 
to the z axis ("axial superlayers"). The remaining four superlayers of six wires each are 
tilted 3° in the direction ("stereo superlayers"). The combination of axial and stereo 
measurements yields a three-dimensional track. Where appropriate, this track is augmented 
with SVX measurements to obtain precise impact parameters. The momentum resolution 
of such tracks, often simply called "SVX tracks," is (j{pt)/pt = [(0.0009pt)^ + (0.0066)^)] V2 
with Pt in units of GeV/c. With such momentum and impact parameter resolutions, along 
with the narrow beam, CDF at the Tevatron is an excellent tool for the study of B physics. 
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The calorimetry systems of CDF lie outside the tracking systems and solenoid. We 
focus on the calorimetry in the |?7| < 1 ("central") region, which is segmented into rj- 
projective towers covering 15° in azimuth and 0.11 units in rj. The inner layers of the 
towers, which make up the central electromagnetic calorimeter (CEM) consist of a 
lead-scintillator stack 18 radiation lengths deep. The CEM has a resolution of (j{Et)/ Et = 
[(0.137)VEt+ (0.02)2]i/2^ where Ft = E sin 9 and E is the measured energy of the tower in 
GeV. A layer of proportional strip chambers (CES) is embedded in the CEM near shower 



maximum and provides measurements of shower position and profile in azimuth and z []17 



The outer layers of the calorimeter tower, which make up the central hadron calorimeter 
(CHA), consist of an iron-scintillator stack 4.5 interaction lengths deep and yield an energy 
resolution of a{ET)/ET = [(0.50)7^^ + (0.03)2]V2. in this analysis, the CHA is used 
primarily to distinguish electrons and photons, which are typically absorbed in the CEM, 
from hadrons, which typically deposit most of their energy in the CHA. 

A three-level trigger system is employed at CDF to select pp events of interest The 
first-level trigger relevant to this analysis selects events based on energy depositions in logical 
"trigger towers" which consist of two adjacent (in 77) calorimeter towers. The second-level 
trigger forms clusters of trigger towers. This trigger level also incorporates a hardware track 
processor (CFT, Central Fast Tracker) |19|, which searches for tracks in the CTC using hits 
in the axial layers and matches those tracks to calorimeter clusters. The third-level trigger 
uses software based on optimized offline reconstruction code to analyze the whole event. 
Details of the trigger selection are given in the next section. 



III. DATA 

The data used in this analysis were collected with triggers which selected events with 
calorimeter signatures characteristic of electrons and photons. During most of the 1994- 
95 data-taking period ("Run IB"), the first-level trigger selected CEM trigger towers with 
minimum Et of 8 GeV. The cross section of this trigger was ~ 20 /xb. 
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Subsequent filtering of tlie surviving events was performed with the specialised "penguin 
trigger", which is a collection of requirements on all three final products of the penguin 
decay chains 5^ —>■ K*'^{^ K~7r^)'~f and 5^ —>■ <!){—>■ K~^K~)'j. The innovative feature of 
this trigger was the use of all the information available at the second trigger level to select 
a specific topological configuration of the final state particles. 

The second-level trigger performed tower clustering and required the event to contain a 
cluster with Et > 10 GeV in the electromagnetic section. The same cluster could include 
hadronic energy deposition and the trigger required the hadronic component to be less 
than 12.5% of the electromagnetic component. A further requirement of at least 4.5 GeV 
deposition in the CES reduced the trigger rate by half while keeping 90% of the electrons 
and photons. 

The CFT track processor was then used to select topologies suggestive of a penguin 
decay, with its photon and two charged hadrons. No track found by the CFT was allowed to 
point at the the same (p as the photon calorimeter tower (spanning 15° in 0). Two oppositely 
charged tracks with px > 2 GeV/c were sought close to the photon (within two calorimeter 
towers) and they were required to lie within 18° of one another in (p. Figure illustrates 
the trigger topology. These track-related requirements were ~ 35% efficient for selecting 
penguin events while reducing the trigger cross section to ~ 80 nb. 

When the trigger rate exceeded the limit of the data taking rate we further reduced the 
trigger rate by rejecting some fraction of the events which satisfied the trigger requirement 
( "prescale" ) . The second-level trigger was prescaled by a factor of two whenever the in- 
stantaneous luminosity was above ~ 21 x 10^° cm^^s^^. The data loss due to the prescale, 
however, was minimal: this trigger considered (22.3 ± 0.9) pb^^ out of the ~ 23 pb^^ of 
data available to it. 

Events satisfying the second-level trigger were then passed to the third-level trigger 
for further consideration. The photon candidate's electromagnetic Et, reevaluated with 
clustering software, was required to be at least 7 GeV, with an associated hadronic energy 
deposition of no more than 15% of that in the GEM. The profiles of energy deposition in 
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the CEM and CES were also required to be consistent with expectations based on test beam 
results for electrons. The track cuts applied by the second-level trigger were confirmed at 
this trigger level using offline beam-constrained tracking in the CTC 

The open points of Figure |^ show the penguin trigger rates as a function of instantaneous 
luminosity during Run IB. These rates can be compared with the total trigger rates at each 
trigger level, shown by the closed points. From this figure we see that one out of 200 events 
accepted by the generic level-one calorimeter trigger also satisfied the second-level penguin 
trigger. The third-level trigger requirements provided an additional rate reduction by a 
factor of 6.5. Approximately 300000 events were collected during Run IB by the penguin 
trigger. The overall trigger efficiency for penguin decays resulting from B mesons with 
Pt > 12 GeV/c and \y\ < 1.25 was (1.7 ± 0.2)% for Ba K*^-i and (2.6 ± 0.3)% for 
Bs 07 decays. In the data sample collected by the penguin trigger in RunIB we expect 
around 7 Bd ^ K*^'y and 2.6 Bg 07 events. This sample was further refined in the offline 
analysis by selecting photon candidates in the good fiducial areas of the calorimeter, and by 
requiring that full CTC track reconstruction revealed no three-dimensional track pointing 
to the cluster. The -E't(7) threshold was raised to 10 GeV. The hadronic/electromagnetic 
energy ratio selection was tightened to 10%, and requirements on shower profile consistency 
were also tightened. 

The trigger thresholds for the penguin trigger were lowered for the 1995-96 data-taking 
period ("Run IC"). At the first trigger level, the Ex threshold was lowered to 5 GeV, 
raising the cross section to ~ 30 fih. The second-level energy requirements were lowered 
to 6 GeV in the CEM and 3 GeV in the CES while the relative hadronic energy and track 
topology requirements were kept the same. The trigger cross section at this level was thus 
raised to ~ 500 nb. The photon Ej- threshold was lowered to 5 GeV in the third- level 
trigger, while the other requirements were kept the same as in RunlB. Due to the lower 
photon energy requirements, the Run IC trigger acceptance rate was six times higher than 
the RunIB trigger, and the signal yield increased by a factor of five. As a result of these 
adjustments, approximately 500000 events were collected from the only (6.6 ± 0.3) pb~^ of 
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Run IC integrated luminosity. The offline Et cut was accordingly lowered for this data to 
8 GeV. 

A sample of electron candidates was also accumulated through Runs IB and IC. The 
trigger for this sample used the same first- level requirements as described above, but required 
Et > S GeV at the second level, along with a CFT track with pr > 7.5 GeV/c pointing to the 
EM cluster's bin. At the third trigger level, the reevaluated thresholds were Et > 7.5 GeV 
and pt > 6 GeV/c. Moreover, the track's trajectory was extrapolated to the CES and 
compared with the shower positions; agreements within ±3 cm in the azimuthal direction 
and ±10 cm in z were required. These trigger requirements were applied throughout Runs 
IB and IC. 

The electron candidate sample serves two purposes in this analysis. In Method I we 
search for radiative decays among events selected by the penguin trigger. The electron 
sample provides a reference signal, B — > e~D^{-^ K~7r'^)X, which we compare to the yield 
of radiative decay candidates. To facilitate this comparison, the same fiducial, Et, and 
calorimeter requirements were applied offline to the subsample of the electron data which 
was collected concurrently with the penguin trigger; the uncertainties in the integrated 
luminosities of these two data sets are thus completely correlated. Because this reference 
sample was obtained by triggering on electrons, a single track was required to point to the 
electron cluster. Nevertheless, in order to simulate the penguin trigger requirements, no 
other track was allowed to point to that bin. 

In Method II, where the photons are identified through their conversion to e+e^ pairs, the 
search for radiative decays is performed in the electron candidate sample itself. In this case, 
the offline selection applies fiducial, shower profile, and track-shower match requirements 
in a manner similar to Method I, but the Et threshold is lower at 8 GeV. The minimum 
track pt is 6 GeV/c. The hadronic/electromagnetic energy ratio requirement is tightened 
to 4% when only one track pointed to the cluster, but is left at 10% in cases with more 
than one track associated with the cluster. This sample also provided the reference signal, 
S+ J/ip{-^ e'^e~)K~^, and thus the entire Run IB data set is used for this method. 
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The electron trigger accumulated 74 pb during this period, amounting to approximately 
3 million events satisfying the offline criteria. 

IV. METHOD I: PHOTON TRIGGER 

In this section, we describe the search for B^^ K*^{-^ K~7T^)'~f and 5^ —>■ 
K~^K~)'~f decays using the penguin trigger described in the previous section. The sensitivity 
of this method to A° ^{~^ pir^)'-/ is strongly reduced by the trigger requirement of 
Pt > 2 GeV/c for the pion track, because in the A — * pvr decays the proton carries most 
of the momentum of its parent and the pion is very slow. Thus, we do not attempt to 
reconstruct such decays. We derive the branching fraction limits for the radiative B decays 
from the ratios between the numbers of candidate events and events of the reference signal, 
B D^{-^ K^Ti^)X, found in the single electron data set. 

A. Radiative Decay Reconstruction 

— *o 

We selected candidate daughters of the K and (f) mesons from the radiative B decays 
by asking for two oppositely charged tracks reconstructed with the inclusion of at least 
three hits in the SVX. Each track was required to have been found by the trigger system 
and have > 2 GeV/c. The penguin trigger topology requirements on the tracks and the 
photon candidate were reinforced offline. We then constrained each pair of candidate tracks 
to intersect at a common vertex and required the confidence level (CL) of the constrained 
fit to exceed 1%. 

We retained two-track combinations consistent with K K^n'^ by requiring 

|M(K^7r+) — M-^.o| < 80 MeV/c^, where M^*o is the world average K*^ mass 
(896.1 MeV/c^) ||^. This window, corresponding to three times the natural K*^ width, con- 
tained more than 85% of the K signal. If the track pair also fell within the mass window 
when the K and vr mass assignments were switched, we chose the assignment which yielded 
the two-track mass closer to the world average. This approach yielded the correct assignment 
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88% of the time. For -> K+K' decays, we required \M{K+K-) - M^\ < 10 MeV/c^ 
where is the world average (p mass (1019.4 MeV/c^) This window, corresponding to 
four times the natural (p width, contained 86.5% of the (p K^K~ signal. 

In order to reject -k^-k^ decays, we assigned pion masses to the two tracks and 

required that |M(7r+7r^) — M^o| > 15 MeV/c^. We thus rejected combinations with masses 
within 2(7 of the world average mass and retained 95.4% of the K*^ K^tt^ decays 
and all of the K^K~ decays. 

The track pair was combined with the photon candidate by adding their four-momenta. 

The trajectory of the photon candidate was determined by assuming that it originated from 

the pp vertex closest in z to the track pair vertex; we call this pp vertex "primary." Because 

— *o 

the lifetimes of the and K mesons are almost ten orders of magnitude smaller than 
that of the B meson [Q, the common fitted vertex of the two charged tracks indicated the 
point where the parent B meson decayed. We computed the B meson's signed decay length 
Lt = Vt ■ Pt/pti where Vt is the displacement in the transverse plane of the B decay 
vertex with respect to the primary vertex (see Figure^), and Pt is the B meson momentum 
projected on the same plane. The proper decay length ct could then be calculated with 
ct = Lt ■ (M/pt), where M is the reconstructed mass of the B meson candidate. The 
typical ct resolution was 100 fim. We required < ct < 3 mm, which retained 90% of the 
signal while rejecting half of the fake B meson candidates formed by tracks coming directly 
from the primary vertex. 

We further required that the B meson carry most of the momentum in its vicinity. We 
defined the isolation variable 

Ib ^ P^(^) (1) 

^ Pt{B) + Ear<iPt' 

where the sum is over tracks consistent with originating from the primary vertex and within 



AR = y {Arjy + (A^)^ < 1 of the B candidate trajectory. The B candidate daughters were 
excluded from the sum. We required Ib > 0.7. Studies with reconstructed B decays in data 
indicate that this requirement is ~ 95% efficient in selecting real B mesons of px > 15 GeV/c 
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while rejecting half of the combinatorial background. 

The mass resolution of B mesons reconstructed in the above manner is given by sim- 
ulation to be 110 MeV/c^, dominated by the energy resolution of the photon. We have 



simulation closely reproduces the momentum resolution and impact parameter resolutions 
of tracks, as well as the energy resolution and shower characteristics of electromagnetic ob- 
jects. After the above selection criteria, there are ~ 400 K 7 and ~ 40 07 events within 
220 MeV/c^ of the world average 5^ and masses of 5279 MeV/c^ and 5369 MeV/c^ 
respectively 0]. 

To further improve our sensitivity to the radiative decays, we exploited the long B meson 
lifetime and the fact that we reconstructed all its daughters. The long lifetime resulted in 
large impact parameters for the K and daughters with respect to the primary vertex; we 
cut on the significance of the impact parameters in the transverse plane, \d/(Jd\. The impact 
parameter resolution was typically ~ 30 fim. We also formed an "alignment angle" 
between the transverse momentum pt and the displacement Vr of the B meson candidate 
(see Figure ^): 



Since we fully reconstructed the B meson, real mesons yielded small values of t^aiign, whereas 
the combinatorial background peaked away from zero. As a pure background sample we used 
events in the high mass region 6 < M(K 7,^7) < 10 GeV/c , where no real B mesons 
should be found. Comparing the OaUgn distributions of the simulated signal events with the 
distribution obtained from the background sample, we selected signal-like events by demand- 
ing Oaiign < 0.15 rad, for both the B^ and B^ decays. We subsequently found the impact 
parameter significance cut which gave the highest signal-to-background efficiency ratio. It 
turned out that the best value was the one which rejected all events in the background (high 
mass) region. 

The optimized selection cuts for 5^ radiative decays were t^aiign < 0.15 rad and 



used D' 



K 7r+ and electrons from the reference signal B ^ e D^X to verify that the 




(2) 
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\d/(Jd\ > 5. These requirements were 66% efficient in retaining K -y decays. 

For the decays, the narrower resonance, compared to the K , resulted in a smaller 
number of combinatorial background events falling within the, consequently narrower, mass 
window used to select the relevant two-track pairs. Thus, the optimized \d/ad\ cut for the 
was less strict at \d/ad\ > 2.5. These optimized requirements are 69% efficient in retaining 
-^s ~^ 4'1 decays. 

All together, the offiine analysis requirements were 14.7% (8.3%) efficient in selecting 
B^ K ■J decays which already satisfied the trigger requirements in Run IB (Run IC). 
The respective efficiencies for 07 decays were 12.9% in Run IB and 7.2% in Run IC. 

The Run IC efficiencies were lower than Run IB because the offiine cut for the photon energy 
was placed at Et > 8 GeV, instead of the 6 GeV trigger threshold, in order to match the 
energy threshold of the electrons used to reconstruct the reference channel. 

Figure |] shows the invariant mass distributions of the three-body combinations surviving 
all the selection criteria. The ±220 MeV/c^ signal region around the world average B 
mass is double-hatched in the figure, and the sideband regions, 3.9 < M{K 7) < 4.9 and 
5.7 < M{K 7) < 6.7 GeV/c^, are single-hatched. One B^ K ^ candidate, from the Run 
IC sample, remains in the signal region, while five populate the sidebands. The expected 
background in the signal region, assuming a uniform distribution interpolated between the 
sidebands, is A^bg = 1-1 ± 0.5 events. There are 2 events just outside the signal window. 
However, the probability of them being signal is small. 

In the 5^ 07 case, no candidates survive the selection cuts. Since there are also no 
events in the i?^ sidebands, in the signal region we expect A^bg < 0.54 events with 90% 
confidence [Q, assuming a uniform distribution interpolated between the sidebands. 

B. Reference Signal Reconstruction 

We reconstructed our reference sample of S ^ e~D^{-^ K~tx^)X decays, by adding the 
four- momenta of the two tracks and the electron candidate. For e + D^ combinations from B 
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decays, we expected the kaon from the to have the same charge as the electron. The mass 
assignment of the pion and kaon masses to the two tracks was thus uniquely determined. 

We retained B D^{-^ K^'k^)X candidates with a Pt{B) distribution similar to that 
of the radiative decay candidates by requiring pTieKn) > 15 GeV/c in Run IB. For Run IC, 
this threshold was lowered to 13.5 GeV/c to accommodate the lower photon threshold. We 
also required that the mass of the three-body combination M{eKTT) be less than 5 GeV/c^. 
Finally, we applied the same < ct < 3 mm and Ib > 0.7 requirements as on the radiative 
decay candidates. These semileptonic decays, however, were not fully reconstructed, and 
we used the combined momentum of the e + system for the (pseudo-proper) lifetime 
calculation. In addition, rather than extrapolating the decay vertex to the trigger elec- 
tron track in order to locate the B decay vertex, we simply used the decay vertex for 
the calculation of ct to avoid additional systematic uncertainties due to the further vertex 
reconstruction. 

We then required {d/cxdl > 3 for the kaon and pion tracks from the —>■ K~t[^ decay. 
Since B D^X decays are not fully reconstructed, we do not make a "i^aiign cut. The 

invariant masses of the selected K~tt^ combinations from B —>■ e~ D^X candidates are shown 
in Figure |^. The K^t^^ combinations with the wrong charge correlation with the electron 
are also shown. We estimated the number oi B e^D^X candidates by fitting the data 
with a Gaussian signal and a linear background and we found 40.7 ± 7.3 events in Run IB 
and 27.4 ± 6.2 events in Run IC. 



C. Efficiencies 



In method I we infer the radiative decay branching fraction from a measurement of its 

ratio with the known B{B —>■ e^D^X). The 6-quark production cross section cancels in the 

— — *o 

ratio, while the effect of systematic uncertainties is reduced. We write for 5^ ^ K 7 , 

B(bI ^ K*%) = B(B ^ e-D^X) x 



DO 
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(3) 



and for 



X 



B(bI (f)-f) = B{B e-D^X) x 

^(j>'y -L/peng 



<t>i 



NeDO 
1 



(4) 



Ju + fdBiD^^K-Tr+)e,DO L^x 
where N^'O^^^^/ N^do is the ratio of the observed number of events of the radiative decays 
and B e^D'^X, eK*o-y,(j,-y/ Cg^io is the ratio of the efficiencies, and i^peng/ Lex is the ratio 
of the integrated luminosities of the penguin and the inclusive electron data samples. We 
assume that the composition of i? — e^D^X candidates is only B~ and 5^, and thus the 
ratios of the fragmentation fractions are fd,s/{fu + fd), neglecting the small contributions 
from other b hadrons such as 5^ and to the denominator. We note that the contribution 



of the 5° through the 5^ — >• e i^e-D** e D^X decay is estimated to be less than 3% in 



the e + D^ sample. The branching fractions ^ and fragmentation fractions |20| used in this 
analysis are listed in Table |. 

Since we use electron trigger data collected concurrently with the penguin trigger data, 
the integrated luminosities of pp collisions are the same for the two data sets. The effective 
integrated luminosities of each data set, however, are different for the two due to the different 
prescale factors. The true integrated luminosities for the penguin and electron data set are 
22.3 pb~^ and 16.2 pb^^, respectively, in Run IB, and 6.6 pb^^ and 4.2 pb~^ in IC. We 
assume that all the uncertainties cancel in the ratio. 

The efficiency ratios were evaluated using a combination of simulation and data. We em- 
ployed a Monte Carlo simulation of events with a single b quark to calculate the efficiencies 
of the kinematic and topological requirements imposed on the data. In this simulation, the 
b quarks were generated with a rapidity and momentum distribution based on a next-to- 
leading order QCD calculation that used the MRSDO parton distribution functions ||22 
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and a renormalization scale of /i = /io = \J'm'l + p'^, where rub = 4.75 GeV/c^ is the mass 
of the b quark and pt is its transverse momentum. These b quarks were subsequently 
hadronized into B mesons using the Peterson fragmentation function |^ with a fragmenta- 
tion parameter et = 0.006. The resulting B mesons were then decayed through the channel 
of interest using the QQ Monte Carlo program to model the phase space, helicity, and 
angular distributions of the decay products. 

For the reference channel, we generated different samples for each of the contributing 
decay chains: B e'VeD^; B e-VeD*{^ D°X); B e-VeD**{^ D^X); and B 
e'Ve{DmT)nr followed by {Dmi)^^ D^X, where (i5n7r)nr indicates a D meson produced 
in non-resonant association with extra pions. We then mixed these semileptonic samples 
according to their relative abundances and selection efficiencies to create a representative 
B e~ D^X sample. We fed these events through the detector and trigger simulations 
to obtain the efficiencies. We also used this simulation to calculate the relative effects 
of the photon/electron trigger cuts, the offiine quality cuts, and the track reconstruction 
in the SVX. We considered simulated SVX track reconstruction since the SVX simulation 
incorporated the same hit efficiencies and pattern recognition as the data. 

Second-level trigger efficiencies were studied using data. The efficiency of the CES energy 
requirement was parameterized as a function of electron or photon Ex by analyzing electrons 
in a very pure sample derived from photon conversions. Applying this parameterization to 
the Monte Carlo samples, we find all the efficiencies to be around 95%. The efficiency ratios 
are therefore near unity, and the 2% uncertainty in the ratio is included in the systematic 
uncertainty. 

The efficiency of the CFT trigger requirements for kaons and pions was determined as 
a function of track pt- We found the CFT is 50% efficient at 1.9 GeV/c and 90% efficient 
at 2.4 GeV/c. The efficiency function of the CFT trigger requirements for the electron in 
the reference signal was determined using a heavily prescaled electron data set with a lower 
energy threshold and no CFT requirement; 50% efficiency is reached at 6.0 GeV/c and 90% 
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at 10.0 GeV/c. The plateau efficiency is 0.915 ± 0.010. These efficiency parameterizations 
were apphed to the Monte Carlo samples to study the effect on the ratios of efficiencies. 
The offiine CTC tracking efficiencies for kaons and pions were estimated by embedding 



Monte Carlo-generated tracks into real J/ip ^ events The efficiency rises with 

Pt in the range 200 < px < 400 MeV/c, and plateaus at a value which depends on the 
instantaneous luminosity and the charge of the track. The integrated efficiency for tracks 
with pt > 400 MeV/c is 0.96 ± 0.02. Again, we applied the efficiency parameterization 
to Monte Carlo samples of the decays of interest. For K K~t:^ and — K^K~ 
decays with the requirement px > 2 GeV/c for the kaons and pions, the efficiency of offiine 
CTC tracking was found to be 0.94 ± 0.04. The corresponding efficiency for the K~tt~^ 
combinations from the D° decays is ~ 1% lower due to the lower pt of the tracks. The 
uncertainties in these efficiencies are dominated by the instantaneous luminosity dependence 
of the tracking efficiency and thus cancel in the efficiency ratio. The offiine tracking efficiency 
for the trigger electron in the reference signal was estimated using an independent electron 
data sample to be 0.99 ± 0.01. We therefore estimate the ratio of tracking efficiencies for 
both if*°7 and 07, relative to the reference signal, to be 1.02 ± 0.02. 

The effect of the isolation requirements for the trigger photon or electron, as well as the 
Ib > 0.7 cut for the B meson, depends strongly on the environment of the B decay (e.g., 
b fragmentation, or multiple pp interactions). We expect similar environments around the 
B mesons in the reference and radiative decay processes and consequently the efficiencies 
are nearly equal. Small differences can be expected due to the extra particles produced in 
B e~D^X decays and because the reference signal contains B~ mesons along with B^. We 
simulated the full pp bb environment using the PYTHIA Monte Carlo generator, tuned to 
match the underlying charged particle distributions in i? i^D^X data p6[. We fed these 



events through the detector and trigger simulations and found that the isolation efficiencies 
are somewhat higher for the radiative decay channels than for the reference signal; the ratio 
is 1.04 ± 0.02 for B^ ^ and 1.06 ± 0.02 for 5° ^ (pj. 

Taking all the efficiencies into account, we find that the efficiency ratios between the 
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radiative decays and the reference channel are eK'O^/eeDO = 2.7 in Run IB and 2.0 in Run 
IC. In the 07 case, we find these ratios to be 3.5 in Run IB and 2.5 in Run IC. 

Table | summarizes the elements of the branching fraction calculation for each of the 
decay modes investigated here. The table also shows the "single event sensitivity" S for the 
two penguin decay modes. S is defined here as 

^ branching fraction 

observed number of events 

and can be rewritten with the known quantities by using Eqs. ^ and This quantity rep- 
resents the branching fraction which would result in an average of one event being observed 
in this analysis. The difference in the sensitivities between the — * K*°7 and 5° — > 07 
decay modes is dominated by the difference of the h quark hadronization fractions. 



D. Systematic Uncertainties 

Table |I| lists the sources of systematic uncertainty considered in this analysis. The 
largest contribution to the total is the uncertainty on the yield of i? ^ D^{-^ K^7i^)X 
decays, which is 18% in Run IB and 23% in Run IC. The second largest contribution arises 
from the uncertainty in the measurement of fs/ {fu + fd) f^^l) followed by the uncertainty in 
the product of branching fractions B{B e^D^X) ■ B{D^ — > K^-k^). 

The last significant contribution to the systematic uncertainty comes from the fraction of 
the time when the meson from a B decay is not an immediate daughter of the B meson but 
is instead a decay product of an intermediate excited D state. Depending on how far down 
the decay chain of the B meson the appears, the kinematics of the resulting kaon and 
pion, and hence the reconstruction efficiencies, are different. In the Monte Carlo simulation 
used to determine the efficiency ratios, the nominal fractions of mesons coming from D** 
mesons and (Dmr)^,. states (/**), from D* mesons (/*), and directly from the B meson (/) 
were f** : f* : f = 0.35 : 0.53 : 0.12 ||]. These fractions were varied to 0.24 : 0.62 : 0.14 and 
0.47 : 0.43 : 0.10. We observed a 12% variation in the efficiency in Run IB and 11% in Run 
IC. We take these variations as the systematic uncertainties in the efficiency ratios. 

22 



The rest of the systematic uncertainty contributions have httle effect on the total, which 
is about 30%. For instance, the Monte Carlo efficiency estimates depend on their input 
distributions, such as the px distribution of the incident particles. We re-weight the Monte 
Carlo Pt{B) distribution which is used as the simulation input by the ratio of the measured 
B production cross section to the theoretical prediction. Even though the efficiencies 
for individual channels vary by as much as ~ 20%, the ratios of efficiencies do not change 
by more than 5%. 

Another relatively small effect is the uncertainty in the difference in trigger efficiencies 
for photons and electrons. The difference resulting from the different Et spectra of the 
photons and electrons is accounted for in the Monte Carlo calculation; moreover, we confirm 
that the detector simulation indeed reproduces the characteristics of the electromagnetic 
shower profile using B e^D^X decays in data. We nevertheless assign an uncertainty 
due to the Et differences between the reference channel electron and the radiative decay 
photon to allow for uncertainties in the simulation of the electromagnetic energy clustering 
at the trigger level. We study the effect of varying the relative efficiency by re-weighting the 
photon and electron Ex distribution in the lowest 10 GeV, away from the efficiency plateau, 
by as much as a factor of two (e.g., the weight is applied for 10 < Erp < 20 GeV in Run 
IB). No weighting is applied for energies in the plateau region. Such a modification of the 
threshold induces a change in the individual event rates by as much as 50%, but the ratio 
varies by only ~ 8%, which we take as the systematic uncertainty. 

The efficiency of the CES trigger requirement itself is measured with an uncertainty of 
~ 1.5%. Assuming that the efficiency for electrons is uncorrelated with that of the radiative 
decay photons, we obtain a conservative 2% systematic uncertainty from this source. 

The CFT efficiency was measured with an uncertainty of ~ 1.5% for kaons and pions, 
and 1% for electrons. Due to the spatial proximity of the two tracks in the radiative decays, 
we consider their efficiencies to be 100% correlated and thus assign a 3% uncertainty for the 
efficiency ratio. Another 2% uncertainty comes from the CTC tracking efficiency, 2% from 
the differences in the isolation efficiencies, and 2% from the finite size of the Monte Carlo 
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samples used to calculate the efficiency ratios. 

The uncertainties listed above were combined in quadrature to obtain the total systematic 
uncertainties on the branching fractions of the radiative decays. As shown in Table |T|, the 
total is ~ 30% for and slightly higher for 5°. 

We combine the Run IB and IC systematic uncertainties by assuming that the uncertain- 
ties due to the statistics of the e + D^ candidates and Monte Carlo samples are uncorrelated, 
and any other sources are fully correlated. The uncorrelated systematic uncertainties are 
added in quadrature, while the fully correlated ones are simply added. The total systematic 
uncertainties are 25% for 'B^^ and 31% for 5^ radiative decays. 



E. Results 

Since we observe no significant signal for either 5^] or 5° radiative decays, we set upper 
limits for their branching fractions. We use a conservative procedure which ignores possible 
background contributions to the observed event yields. 

First, we calculate an upper limit on the mean number of radiative decays Nql at a given 
CL, including the total systematic uncertainty CTgyst, by numerically solving the following 
equation: 

1-CL=Y.Vm^,,.^^M, (6) 

n=0 

where A^'obs is the number of candidates observed, and 'P^,o-(^) is defined with the Poisson 
distribution Pfj,{n) and the Gaussian distribution G^^o-(x) as follows: 

POO 

P^,,(n) = / Pxin)G^A^)dx. (7) 
Jo 

With one K 'j candidate observed in the entire data sample and a 25% uncer- 

tainty, the upper limit on the mean number of radiative decays is 4.3 (5.5) at 90% (95%) 
CL. This result, with a single event sensitivity (Equation^) of 3.8 x 10~^, yields upper limits 
on the branching fraction B(bI ¥*%) of 1.6 x 10"^ at 90% CL and 2.1 x 10"^ at 95% 
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CL. With no — > 07 candidates and a total uncertainty of 31%, we expect less than 2.6 
(3.6) events on average at 90% (95%) CL. With a single event sensitivity of 9.3 x 10^^, we 
thus obtain 13(^, 07) < 2.5 x 10"^ at 90% CL and < 3.3 x lO'^ at 95% CL. 

V. METHOD II: PHOTON CONVERSION 

In this section, we describe the search for B^—^K {-^ K^n^)^, B^ 0(~*' K~^K~)'y, 
and A° A(^ pTT~)j decays in which the photon is identified by an electron-positron 
pair produced through photon conversion before reaching the CTC volume. A conversion 
daughter with Ex > 8 GeV served as the trigger; the same inclusive electron trigger was 
used for the e + sample in Method I. 

Though the typical photon conversion probabihty was 6% for CDF in this data, this 
analysis benefits from the fact that we can utilize all of the Run IB data, which corresponds 
to an integrated luminosity of 74 pb~^, or three times more than that collected with the 
penguin trigger, and that there was no requirement of any additional tracks at the trigger 
level. This fact allowed us to apply, in the offline selection, apr threshold as low as 0.5 GeV/c 
to the hadron tracks coming from the b hadron decays instead of the 2 GeV/c cut used in 
Method L This lower threshold essentially doubles the efficiency for the B hadron decay 
products. Moreover, in the relatively low energy region of our interest where the tracking has 
better resolution than the calorimetry, reconstructing b hadron masses from the momenta 
measured by the tracking detectors has the advantage of good mass resolution. This is 
typically 45 MeV/c^ for the reconstructed B mesons and is dominated by the momentum 
resolution of the trigger electron. 

We derive the branching fractions for the radiative b hadron decays from the ratios 
between the numbers of such decays and B~ J/il>{-^ e'^e~)K~ decays found in the same 
data set. The uncertainties in the b quark production cross section and on the integrated 
luminosity thus cancel, as well as most of the uncertainties on the detection efficiency. It 
would have been preferable to use 5^ J/ipK*^, B^ — > J/ip(f), and A° — > J/ipA decays 



25 



instead of B~ J/ipK', since they arise from the same production mechanisms as the 
corresponding radiative decays and are topologically more similar. However our samples of 
those final states are too small to be useful as normalization. 



A. Radiative Decay Reconstruction 

Reconstruction of the radiative decays began with identification of a photon conversion. 
A photon conversion candidate was formed by the electron candidate and an oppositely 
charged track with > 0.5 GeV/c. A fit was made which constrains the two tracks to 
originate from a common vertex and be parallel to each other at the vertex. The CL of the 
fit was required to be greater than 0.1%. The background due to misidentified electrons and 
combinatorial backgrounds is small (< 1%) among the photon conversion candidates with a 
vertex outside the beam pipe. The candidates that have their conversion points inside the 
beam pipe are dominated by real electron-positron pairs from Dalitz and rj decays. We 
required the transverse distance of the conversion point from the nominal beamline to be 
less than 30 cm in order to ensure that it is in the well known materials before the CTC, 
and to be greater than 3 cm in order to reject backgrounds from Dalitz decays. We obtained 
~ 850000 photon conversion candidates in the Run IB data. Figures ^ and || show, for all 
transverse distances, the reconstructed conversion vertex density in the x — y plane and r — z 
plane. The fine structure of the CDF tracking detectors such as the SVX (r ~ 5 cm), the 
VTX (r ~ 15 cm), and the CTC (r > 30 cm) can be clearly resolved. The detailed study 
of the CDF material distribution using ~ 200000 conversion candidates in 1992-1993 data 
is described in ||2^ . 



Q 

For each photon conversion candidate in an event, we searched for ^ K ■y and 
-^s ~^ 4'1 decays. A !b|] candidate was formed by the photon conversion candidate and a 
pair of oppositely charged tracks. The two "meson tracks" were required to be reconstructed 
in the SVX with hits in at least 3 layers. In addition, the transverse momenta had to exceed 
0.5 GeV/c for each track and 2 GeV/c for the two-track system. A fit was performed with 
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the following topological constraints: (1) the meson tracks originate from a common vertex; 
(2) the photon conversion candidate points back to the meson decay vertex; and (3) the four- 
track system points back to the primary vertex, which was defined to be the pp collision point 
nearest in z to the trigger electron track's closest approach to the beamline. We required 
the CL of the fit to be greater than 0.1%. The 5° candidate was then accepted if the 
reconstructed K mass was within ±80 MeV/c^ of the world average value. Both K^tt and 
Tc^K~ mass assignments were considered for the K candidate, and the assignment giving 
a value closer to the world average was chosen. We also required that the pseudorapidity of 
the B candidate \r]B\ be less than 1. Finally, we selected candidates with lifetime ct > 100 
/im and > 0.7 (See Section [IV A| ). 

The selection of -B^ candidate proceeded on similar lines, except both tracks were assigned 
kaon masses and the mass window was ±10 MeV/c^ around the world average. 

At this point, there were 15 i^" 7 and one 07 events within 100 MeV/c of the corre- 
sponding world average B masses. As previously noted, the mass resolution of the recon- 
structed B mesons is about 45 MeV/c^. We refined this selection by tightening the pt cut on 
the two-track system and by applying impact parameter significance cuts to the individual 
meson tracks. The thresholds were optimized by maximizing egig/ y/e^^, where esig and ebg 
are the efficiencies for the signal and background events found in the ±100 MeV/c^ window 
around the B masses. The signal efficiency was obtained from Monte Carlo calculations 
similar to that of Method I (see Section [IV C|) , while ebg was estimated by interpolating the 
observed yields in the mass sidebands, defined to extend from 200 to 1200 MeV/c^ above 
and below the average mass, through the signal region. For the B^^ channel, the optimized 
selection cuts were pxiKn) > 2.75 GeV/c and \d/ad\ > 4.5 for both meson tracks. Fig- 
ure ^(top) shows the K 7 mass distribution after these cuts. Any further cuts, for example 
on the proper decay length, did not improve esig/ y/^- One candidate remained in the -B° 
signal region; the expected background is 0.6 ± 0.3 events. 

For the 5° channel, the optimized selection cuts were pt{KK) > 2.25 GeV/c and 
\d/<7d\ > 3.0. The resulting invariant mass distribution is shown in Figure 0(top). No 
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candidates were found in the signal region, where we expected a background of 0.1 ± 0.1 
events. 

The decay —>■ A{-^ pTr~)j is topologically distinct from the meson decays. Since the 
A has a long lifetime, with ~ 8 cm, it decays outside the SVX fiducial volume ~ 85% of 
the time, and thus only 15% of the A decays are expected to have associated SVX tracks. 
We therefore first reconstructed A's without using SVX information. The higher-pr track 
of the track pair was assumed to be the proton, and was required to have > 1-5 GeV/c 
while the pion had to have pr > 0.4 GeV/c. The energy loss dE/dx for both tracks had to 
be consistent with expectations. A vertex-constrained fit of the track pair was accepted if its 
CL exceeds 0.1%. Photon conversions, a major source of background for A — > pn' decays, 
were rejected here by eliminating those track pairs which could be fit with the conversion 
hypothesis. Finally, the track pair was accepted as a "CTC-A" candidate if the distance of 
the decay vertex from the nominal beamline exceeded 1 cm. 

If both the proton and pion tracks had at least two SVX hits, the vertex-constrained fit 
was redone using the SVX information. Again, the CL of the fit was required to be greater 
than 0.1%. We also required the SVX layer hit pattern to be consistent with the expectation 
from the reconstructed A decay. For example, if the A decay vertex was between the second 
and third of the four SVX layers, we required that the tracks have exactly two hits in the 
outermost layers. About 10% of the "CTC-A" candidates satisfied the above requirements 
and were thus reclassified as "SVX-A" candidates. 

A A° candidate was formed by a photon conversion and a A candidate. From the CTC-A 
candidates, we reconstructed "CTC-A°" candidates with a constraint that both the A and the 
photon point back to the primary vertex. This constraint improved the A° mass resolution 
from 75 MeV/c^, without the constraint, to 50 MeV/c^. For the SVX-A candidates, however, 
only the photon was constrained to point back to the primary vertex, while the A trajectory 
was required only to point backwards to within 2 cm in z of the primary vertex. The typical 
A^ mass for these "SVX-A^" candidates is also 50 MeV/c^. In both cases, we required the CL 
of the constrained fit to exceed 0.1%. We then recalculated the A mass given the constraints 
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and required that it fell within ±3 MeV/c^ of the world average A mass. The typical A mass 
resolutions are 2.5 MeV/c^ for CTC-AO candidates, and 1.5 MeV/c^ for SYX-A^. 

We improved the sample purity by requiring large impact parameters, recalculated after 
the constrained fit, for the proton and pion tracks. In the SVX-A^ case, the impact parameter 
resolution was good enough to require at least 3(7^ inconsistency with the primary vertex. 
In the CTC-A° case, however, we noted that the proton carries most of the momentum of 
its parent and required only \d/ad\ > 0.5 inconsistency. The pion from A° decay is more 
likely to have a large impact parameter, so we required \d/ad\ > 2. Finally, we selected A° 
pseudorapidity |%o| < 1 and isolation Ib > 0.7, as before. After these selection cuts, we 
found 23 CTC-A^ and 2 SVX-A|) candidates in the ±100 MeV/c^ window around the world 
average A° mass. 

The SVX-A^ candidates were further refined by considering the signed impact param- 
eter of the A's. The sign is defined as positive when the crossing point of the A and the 
A° momenta lies in the hemisphere containing the A|J, as should be the case for real A° 
decays. The typical resolution of the signed impact parameter is 40 fim. Following the same 
optimization procedure as before, we find that a cut value of 70 /im maximizes e^ig/ i/cbg- 
No candidates survived this cut, while the expected background is 0.1 ± 0.1 events. 

Since the CTC-A°'s lack the improved impact parameter resolutions of the SVX, we 
reinforced the kinematic requirements by requiring the px of the A to be greater than 
4 GeV/c. Two candidates remained in the signal region, and the expected background 
is 3.3 ± 0.6 events. Combining the CTC and SVX samples, we found two candidates in 
the signal region with an expected background of 3.4 ± 0.6 events. The invariant mass 
distribution is shown in Figure (top). 

B. Reference Signal Reconstruction 

The reference signal for this analysis method consists of B~ J/ipl^ e^e~)K~ decays. 
K J/ip ^ e~^e~ candidate was formed by the electron candidate and an oppositely charged 
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track with px > 1 GeV/c. We required the partner track to exhibit energy loss in the 
CTC and deposition in the GEM in a manner consistent with being an electron. The two 
tracks were then subject to a vertex-constrained fit, and its CL is required to be greater 
than 0.1%. The dielectron invariant mass distribution is shown in Figure The ratio 
of signal to background S/B is approximately 1/2 in the 2.8 to 3.2 GeV/c^ mass range. 
The backgrounds are mostly combinatorial, involving hadrons misidentified as the partner 
electron. The low-mass tail on the signal is due to photon bremsstrahlung on the electron 
tracks. A fit of the mass distribution with two Gaussians and a second-order polynomial 
yields ~ 8000 J/ip —>■ e+e^ events. 

The J/ip candidates were then combined with a track with px > 2 GeV/c. We required 
that all three tracks incorporate at least 3 SVX hits. We constrained the tracks to a common 
vertex pointing back to the primary vertex and accepted the combination if the GL of this fit 
exceeded 0.1%. We also required that the candidate trajectory fall within the pseudorapid- 
ity range |?7b| < 1, have proper lifetime ct > 100 fim, and isolation Ib > 0.7. The resulting 
M{eeK) mass distribution shows the same low-mass bremsstrahlung tail as the M(ee) dis- 
tribution; in order to correct for it and, at the same time, compensate for the resolution lost 
because of the electron momentum uncertainty, we plot M{eeK) — M{ee) + Mj/^, where 
Mj/^ is the world average J/ip mass, instead of M{eeK). The resolution on this compen- 
sated mass is typically 25 MeV/c^, whereas it is typically 50 MeV/c^ for M{eeK) alone. 

After the above selection, we have 48 J/ip candidates with S/B 10 in the ±100 MeV/c^ 
window around the world average B~ mass. Further requirements, determined by the cut 
optimizations on the different radiative decays, were applied to this sample in order to 
achieve as much cancellation of the systematic uncertainties as possible. To compare with the 
B^^ K 7 decays shown in Figure [^(bottom), these requirements are pt{K) > 2.75 GeV/c 
and \d{K)/ad\ > 4.5. The signal yields were calculated by subtracting the backgrounds 
estimated from the sidebands, which range from 200 to 300 MeV/c^ above and below the 
B^ mass. The yield is 28.0 ± 5.8 events. In the 5° 07 case shown in Figure pUKbottom), 
the cuts are Pt{K) > 2.25 GeV/c and \d{K)/ad\ > 3, yielding 35.0 ± 6.4 events. For 
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the Aj) A7 case, only the Pt{K) > 4 GeV/c cut was apphed. The yield, shown in 
Figure |TT](bottom) , is 24.0 ± 5.3 events. 



C. Efficiencies 



Because no significant excesses over backgrounds were observed in any of the radiative 
decay modes investigated, we set upper limits on the branching fractions. As in Method I, 
we start from the ratios between the number of observed signal and reference decays. Since 
these decays were reconstructed in the same data set, the b quark production cross section 
and the integrated luminosity of the data cancel in this ratio. The fragmentation fractions, 
branching fractions, and total reconstruction efficiencies, on the other hand, do not cancel 
in principle, and their ratios must be estimated. We write the following relations: 
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The branching fractions [^] and fragmentation fractions ||20| which we used are listed in 
Table |T|. The remainder of the calculation concerns the efficiency ratios. The efficiency 
ratios for most kinematic and geometric requirements, including those on E^, Pr, masses, 
ct, impact parameters, and fit constraints, can be reliably calculated with simulation, as 
in Method I. Likewise, the effect of the electron trigger can be calculated by applying an 
efficiency curve as a function of electron Et and pt to the Monte Carlo samples, where the 
curve is based on measurements using unbiased data collected with independent triggers. 



31 



We assume that the B isolation cut efficiencies cancel exactly in the ratio, since, unlike in 
Method I, the reference decay is fully reconstructed. 

The effect of the tracking efficiencies on the ratio is also mostly included in the Monte 
Carlo calculation, but since the radiative b decay leaves four tracks and the reference de- 
cay only three, we accounted for the second meson track by multiplying the Monte Carlo 
efficiency by the integrated CTC tracking efficiency, 0.96 ± 0.02, estimated by embedding 
simulated tracks in CDF data |2^ (see Section [IV Q ). As previously noted, the Monte Carlo 
simulation already models the SVX efficiency, and thus no further correction to the tracking 
efficiency is needed. 

Effects which do not cancel in the ratio include the efficiencies of the quality cuts for the 
J/ip partner electron, the A p'K~ selection, and the conversion probabilities. 

The quality cut efficiency of the J/ip partner electron was estimated from the J/ip can- 
didates themselves to be 0.75 ± 0.03 by counting the number of the J/ip signals before and 
after the quality cut. In a similar manner, the A quality cut efficiency was estimated to be 
0.72 ±0.02. We investigated the effect of the photon conversion probability in detail because 
it dominates the total efficiency differences between the radiative b decays and the reference 
decay. 



The detector simulation, described in Section [IV Q , also simulates photon conversions. 
The material distribution of the CDF inner detector used by the simulation is based on 
previous photon conversion measurements and a careful accounting of the material of the 
CTC inner wall which is known to be (1.26 ± 0.06)% of a radiation length. We calibrated 
the simulation by normalizing the conversions simulated in the CTC inner wall with the 
rate seen in the data. The data used consists of the K 7 candidates, but with loose 

selection cuts on ct, J^, and mass to increase the sample size. The resulting conversion 
probability from the Monte Carlo calculations is ~ 6%. The simulation was analyzed in 
the same manner as the data; in this way, the non-uniformity in the material distribution 
and the consequent dependence of the conversion probability on the physics process and 
event selection criteria was included in the simulation calibration. In particular, requiring 
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the meson tracks to be reconstructed in the SVX, as is the case in the and samples, 

imphes that most of the photons will pass through approximately 1%Xq more material than 

those in events where the tracks lie outside the SVX fiducial volume. On the other hand, 

the A° analysis makes no SVX requirements on the tracks; since 50% of such photons are 

outside the SVX volume, they traverse, on average, ~ 0.5%Xo less material compared to 

the B meson case. The process dependent scale factors which relate the data samples to 

— — 

the simulation normalization are found to be 0.89 ± 0.05 for the B^ and B^ decays, and 
0.95 ± 0.05 for the decay. 



Table |T| shows a summary of the efficiency estimates for each of the decay modes. For 



example, the ratio for B^ K 7 is given by 0.064 x 0.89 x (0.96/0.75), where 0.064 is the 
Monte Carlo efficiency ratio, 0.89 is the conversion probability scale factor, 0.96 is the CTC 
tracking efficiency for the second meson track, and 0.75 is the partner electron quality cut 
efficiency for the J/ip ^ e~^e~ decay in the reference sample. As expected, the efficiency 
ratio is around 6%, largely due to the conversion probability. 



The single event sensitivities defined by Eqs. § and |8|-|T0| are also shown in Table pi 



They are 4.4 x 10 ^ for 5°, 9.5 x 10 ^ for and 2.8 x 10 for A°. The differences among 
the sensitivities are dominated by the differences among the b quark fragmentation fractions. 



D. Systematic uncertainties 



Table summarizes the sources of systematic uncertainties for each of the decay modes 
considered in this analysis. One of the largest uncertainties arises from the statistical un- 
certainty in the J/ipK yield, contributing 21% for 5^], 18% for 'B^, and 22% for the A° 
channel. The uncertainty due to the input branching fractions is dominated by that of 
B{B^ — > J/ipK^), and we assign it 11% for all the decay modes. 

The other major source of systematic uncertainty is the measurement of the fragmen- 
tation fractions fsl fu and fi^ol fu 01- These fractions were measured at CDF using the 
decays 5^ e^DfX and A? e^A'tX, normalized to B^ — > D^X. Their quoted un- 
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certainties are 18% for fs/fu and 35% for /ao//u, but these values include a 6% uncertainty, 
originating from the h hadron pt spectrum, which is fully correlated with the corresponding 
uncertainty in this analysis. We thus reduced the quoted uncertainties by 6% in quadrature 
and obtained a 17% systematic uncertainty due to fs/fu and 34% due to f^o/fu- 

We confirmed that changing the h quark px spectrum does not contribute any systematic 
uncertainty, since this spectrum is common to all the decay modes, by changing the Monte 
Carlo generation parameters from their nominal values = 4.75 GeV/c^ and /i — /iq- The 
b quark mass was changed to 4.5 and 5.0 GeV/c^, and the renormalization scale was changed 
to Hq/2 and Individual efficiencies for the radiative and B~ — > J/ipK' decays vary by 
~ 20%, but the efficiency ratios remain, as expected, stable within the uncertainties of the 
finite Monte Carlo samples. 

Small systematic uncertainties are contributed by efficiency factors which do not cancel 
in the ratio. For instance, for the photon conversion probability correction, which was 
evaluated to be 0.89 ± 0.05 for the B mesons, we assign a 6% systematic uncertainty. For 
the case, the uncertainty is 5%. We assign a 4% systematic uncertainty for the quality 
cut efficiency on the partner electron in the J/ip ^ e^e" decay, and 3% for the quality cut 
efficiency for reconstructing A — > p7r~. These two uncertainties arise from the data sample 
sizes used for the efficiency estimation. The CTC tracking efficiency contributes another 2% 
systematic uncertainty which comes from its instantaneous luminosity and electric charge 
dependence. 

Another effect which does not cancel in the efficiency ratio is that the 
hadronic/electromagnetic energy ratio cut depends on the number of tracks pointing to 
the calorimeter cluster. This number is different for photon conversions and J/V' — > e'^e" 
decays. About 45% of the conversion partners point to the same cluster as the trigger elec- 
tron, while less than 1% of the partner electrons in J/-^ decay exhibit the same behavior. 
In principle, the effect of this difference can be estimated with a full simulation of the pp 
event, including h fragmentation products and multiple pp collisions. Instead, we estimated 
this systematic uncertainty to be about 5% based on the efficiency difference between the 
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two different hadronic/electromagnetic energy ratio cuts on the J/ip ^ e^e~ candidates in 
the data. 

Finally, the systematic uncertainties due to the finite Monte Carlo sample sizes in the 
efficiency calculations were all around 4%. When all these uncertainties were combined in 
quadrature, we found the total systematic uncertainties to be 26% for 29% for 5°, and 
43% for AO. 



E. Results 

The low background level for and 5° radiative decays allows us to set limits on the 
branching fractions without background subtraction. For the case, however, we account 
for the expected background level by using a simple simulation which generates the numbers 
of signal and background events in each trial according to the probability distributions 
'^AfcL,<Tsyst('^) and "PiVbg.abgl"-), where V^^„{n) is defined in Eq. 0. A^cl is the upper limit on 
the number of decays for a given CL, a^jst is the systematic uncertainty on the signal yield, 
and A^bg is the number of background events with uncertainty abg. The CL is given by the 
fraction of trials which has the total number of signal and background events exceeding the 
observed number of events Aobs, but still has fewer background events than Aobs- 

We calculated A'cl to be 4.3 for 2.6 for and 4.5 for A|] at 90% CL, and 5.5, 3.5, 



and 6.8, respectively, at 95% CL. With the single event sensitivities listed in Table 111, we 



obtained the limits on the branching fraction, B(W^ —>■ K*^'-/) < 1.9 x 10 ^ (2.4 x 10 
(95%) CL. 



B{bI ^ 07) < 2.5 X 10-^ (3.4 x 10"^), and B{Al A7) < 1.3 x lO^^ (1.9 x lO^^) at 90% 



VI. COMBINED LIMITS 

Q Q 

Since the two analyses searching for ^ K •y and B^ — > 07 decays are statistically 
independent, we simply add the numbers of candidates found in each analysis. In total, 
there are two 5° candidates with an expected background of 0.6 ± 0.3 events, and no 5^ 
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candidates with an expected background of 0.1 ±0.1 events. The combination does not yield 
any significant excesses over the background level but does tighten the upper limits on the 
branching fractions. 

The combined single event sensitivity of using both methods is given by SY^u — SJ^ + 
Sj'l and is 2.0 x 10~^ for and 4.7 x 10~^ for S^. The systematic uncertainties due 
to the generated Pr(-S) spectrum, fs/fm ^{4' ~^ K~^K~), and CTC pattern recognition 
efficiency are fully correlated between the two methods and simply added together; the 
other systematic uncertainties are considered to be fully uncorrelated and are thus added 
in quadrature. We obtained 18% as the combined systematic uncertainty for 5^ and 25% 
for S^. We then calculated, without any background subtraction, the upper limits on the 
branching fractions B(^a ^ ^*%) < 1-1 x 10"^ (1.4 x 10"^) and B(bI ^7) < 1.2 x 10"^ 
(1.6 x 10-^) at 90% (95%) CL. 

VII. CONCLUSIONS 

We have searched for ^ X*°(^ K-n+)-f, 5^ ^ 0(^ K+K-)-f, A(^ PT^'h, 

and their charge conjugate decays, using events produced in pp collisions at = 1.8 TeV 
and recorded by CDF. Two methods were employed. 

In the first method, the photon was detected in the electromagnetic calorimeter as a 
cluster of energy. We designed and installed a dedicated trigger which, in addition to the 
photons, required information about the charged particles originating from the daughter 
meson. We collected 22.3 pb"^ of data with £^t(7) > 10 GeV during 1995 and 6.6 pb"^ of 
data with Eri'j) > 6 GeV during 1995-96. 

In the second method, the photon was identified by an electron-positron pair produced 
through external photon conversion before the tracking detector volume. One of the con- 
version electrons with Et > 8 GeV served as a trigger for event recording; no additional 
tracks coming from the daughter hadron decay were required. The trigger recorded 74 pb~^ 
of data from the 1994-96 period. We observed no significant signal in both the methods. 
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and set upper limits on the branching fractions (Table 0). 

Combining the two analyses, we obtained upper limits on the branching fractions 



^Z*°7) < 1.4 X 10-^ 
B(bI 07) < 1.6 X 10-^ 
S(A° A7) < 1.9 X 10-=^ 



at 95% CL. The result on the — > i^' 7 decays is consistent with the measurements per- 
formed in the e+e^ colliders The results on the 5^ and decays are the current lowest 
limit on these branching fractions and they are also consistent with the theoretical prediction 
that the 5^ ^7 and 5^ K*^'~f branching fractions are of the same magnitude [lin . 
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FIGURES 



FIG. 1. Feynman diagrams for the b ^ and h ^ d'y penguin loops. 

FIG. 2. The topology of the objects considered by the penguin trigger, shown on a schematic 
depiction of the CEM calorimeter with the beam pipe going perpedicularly through this page. 

FIG. 3. Level 1, 2, and 3 trigger rates for the photon + 2 track trigger as a function of 
instantaneous luminosity in Run IB (open points). Total trigger rates for each stage are also 
shown (filled points). 

FIG. 4. The B decay vertex and relevant quantities on the plane transverse to the beam. For 
clarity, only the B momentum and one of its' charged daughters are shown. 

FIG. 5. Top: 'yK~iT'^ invariant mass distribution for 5° — ^K*^{-^ K~'it'^). There is one 
candidate. Bottom: 'yK~^K~ invariant mass distribution for B^ — > 7^(^ K~^K~). There are no 
candidates seen. 

FIG. 6. Invariant mass distributions of the K~Tr~^ combinations for B — > e~-D°(— > K~Tr~^)X, 
decays in the Run IB (top) and IC (bottom) data. The right-sign distributions (points) are for 
same charge electrons and kaons, as should be the case if they are both products of the real B 
decay chain, whereas in the wrong-sign distributions (histograms) the kaon has opposite charge 
to the electron. By fitting a Gaussian and a straight line to the right-sign distributions we find 
40.7ib7.3 and 27.4ib6.2 candidate B — > e~D^{-^ K^-k^)X events in Runs IB and IC, respectively. 

FIG. 7. Photon conversion vertex density in the x — y plane in the 74 pb^^ of CDF Run IB 
inclusive electron data. The fine structure of the CDF tracking detectors can be clearly resolved. 

FIG. 8. Photon conversion vertex density in the r — z plane in the 74 pb~^ of CDF Run IB 
inclusive electron data. The fine structure of the CDF tracking detectors can be clearly resolved. 
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FIG. 9. Top: e^e'K^'K'^ invariant mass distribution for K (— iC~7r+)7(— e^e~) in 

the 74 pb""*^ of CDF Run IB inclusive electron data. Bottom: corresponding e^e~K~ invariant 
mass distribution for the B~ —>■ J/ijj{—>- e^e~)K~ reference decay. There are 28.0 it 5.8 events 
after background subtraction. 

FIG. 10. Top: e^e~K'^K~ invariant mass distribution for — ^(— K"^ K~)^{-^ e^e") in 
the 74 pb^^ of CDF Run IB inclusive electron data. Bottom: corresponding e^e^ K~ invariant 
mass distribution for the B~ —>■ J/ijj{—>- e^e~)K~ reference decay. There are 35.0 it 6.4 events 
after background subtraction. 

FIG. 11. Top: e^e~p'K~ invariant mass distribution for A° — A(— piT~)^{-^ e+e~) in the 
74 pb~^ of CDF Run IB inclusive electron data. Bottom: corresponding e^e~K~ invariant mass 
distribution for the B~ — > J/'^{^ e^e~)K~ reference decay. There are 24.0 it 5.3 events after 
background subtraction. 

FIG. 12. Dielectron invariant mass distribution of the J/^ — e+e~ candidates in the 74 pb~-^ 
of CDF Run IB inclusive electron data. The number of the J/^ e+e~ events obtained by fitting 
the mass distribution to a function of 2 Gaussians and a polynomial is ~ 8000. 
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TABLES 

— — *o 

TABLE L Ingredients for the calculation of the branching fractions of -B^ K 7 and 
07 decays in Method I according to Eqs. (|3|) and (Q). 







5° ^07 




Run IB Run IC 


Run IB Run IC 


Nohs (events) 


1 





-^bg (events) 


0.9 ±0.4 0.2 ±0.2 


< 0.54 (90% CL) < 0.54 (90% CL) 


N^jjo (events) 


40.7 ±7.3 27.4 ±6.2 


40.7 ±7.3 27.4 ±6.2 


fd,s/ ifu + fd) 


1/2 


0.213 ±0.038 


B(K*^ ^ K-7r+) 


2/3 


— 






n 401 _|_ nno 

U.^iC/X _1_ W.VJwO 


B(B e-D^X) ■ B{D^ - 


-^K-TT+) (2.94 ± 0.4) X 10^3 


(2.94 ±0.4) X 10-^ 


Csignal/ ^eD" 


2.65 2.01 


3.50 2.48 




22.3/16.2 6.6/4.2 


22.3/16.2 6.6/4.2 


Single event sensitivity 


5.94 X 10"^ 10.2 X 10"^ 


1.44 X 10"^ 2.64 X 10-^ 


Combined 


3.75 X 10-^ 


9.29 X 10-5 
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TABLE II. Systematic uncertainties on the branching fractions oi K 7 and — ^7 

decays in Method I. 



B^^K B^^(j)-f 



Source 


Run IB 


Run IC 


Run IB 




Run IC 


eD^ statistics 


18% 


23% 


18% 




23% 


Monte Carlo statistics 


2% 


2% 


2% 




2% 


Composition of e + D° sample 


12% 


11% 


12% 




11% 


Pt{B) distribution 


3% 


3% 


5% 




2% 


CEM Et cut efficiency 


7% 


7% 


8% 




8% 


CFT efficiency 


3% 


3% 


3% 




3% 


CTC pattern recognition 


2% 


2% 


2% 




2% 


XCES efficiency 


2% 


2% 


2% 




2% 


Isolation efficiency 


2% 


2% 


2% 




2% 


fs/ifu + fd) 








18% 




B(B e-D^X) ■ B{D^ K~Tr+) 


14% 






14% 




B{(t> K+K-) 








2% 




Total systematic uncertainty 


27% 


30% 


33% 




36% 



Combined 25% 31% 
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TABLE III. Ingredients for the calculation of the branching fractions of K 7, — > (j)^, 

and A7 decays in Method II according to Eqs. (|8|)-(|lC|). 

^°^^*°7 ^s^'Pl Ag^A7 

iVobs (events) 10 2 

TVbg (events) 0.6 ± 0.3 0.1 ±0.1 3.4 ± 0.6 

Nj/^K (events) 28.0 ± 5.8 35.0 ± 6.4 24.0 ± 5.3 



fd,s,Ai,/ fu 

B{cl) K+K-) 
B{A PTT-) 
B{B- ^ J/i;K- 
B{J/iP 



e ' e 



1 

2/3 



0.426 ± 0.076 



0.491 ± 0.008 



0.236 ± 0.082 



(0.99 ± 0.10) X 10-3 (0.99 ± 0.10) x lO^^ 
(6.02 ± 0.19) X 10-2 (6.02 ± 0.19) x 10"^ 



0.639 ± 0.005 
(0.99 ±0.10) X 10^3 
(6.02 ±0.19) X 10^2 



CTC tracking 
J/ip partner electron 
A quality cut 
Xt{T)ATA)/ Xt(MC) 



0.960 ± 0.020 
0.749 ± 0.028 

0.889 ± 0.052 



0.960 ± 0.020 
0.749 ± 0.028 

0.889 ± 0.052 



0.960 ± 0.020 
0.749 ± 0.028 
0.721 ±0.018 
0.954 ± 0.047 



Csignal/ ^J/ipK 
^signal/ ^J/ipK 



MC 



0.0644 
0.0733 



Single event sensitivity 



4.36 X 10" 



0.0748 
0.0853 



0.0666 
0.0588 



9.54 X 10" 



2.80 X 10"^ 
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TABLE IV. Summary of the systematic uncertainties for Method II. 





J-'d ^ 7 


— 




JlihT^ statistics 


21% 


18% 


22% 






o/o 




l^OIlVcrblOIl prODclDlilljy 


D/O 


o/o 


o/o 


J / pdi Liiei eieciiuii 






4/0 


A aHi/ax 






qO/ 

o/o 


CTC pattern recognition 


2% 


2% 


2% 


HAD/EM 


5% 


5% 


5% 


Fragmentation fractions 


0% 


17% 


34% 


Branching fractions 


11% 


11% 


11% 


Total 


26% 


29% 


43% 



TABLE V. Summary of the branching fraction hmits. 









K 7 


^s^4n 




► A7 


Confidence level 


90% 




95% 


90% 


95% 


90% 


95% 


Method I 


1.6 X 10" 


-4 


2.1 X 10"^ 


2.5 X 10"^ 


3.3 X 10"^ 






Method II 


1.9 X 10" 


-4 


2.4 X 10-^ 


2.5 X 10-^ 


3.4 X 10-^ 


1.3 X 10-^ 


1.9 X 10-3 


Combined 


1.1 X 10- 


-4 


1.4 X 10-^ 


1.2 X 10-^ 


1.6 X 10-^ 


1.3 X 10-3 


1.9 X 10-3 
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pp colllllision 

FIG. 2. The topology of the objects considered by the penguin trigger, shown on a schematic 
depiction of the GEM calorimeter with the beam pipe going perpedicularly through this page. 
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N 



Run IB data 

I I I I I I I I I I I I I I I I I I I I I I I I I I I 

■ 

_ ■ 

10^^ . ■ ■■ ■ ■ (L1 totol) 



cn 



10 ^ 



Q-i □ n ° ° ('-^ colorimeter) 

2 n ° ° 



, • • •• • 



(L2 total) 



i, i,* **★★* ★ * 

(L3 total) 



★ ★ 



o oo o 



o ° 

(S3 



o ^ 



o o (L2 penguin) 



-1 * A (L3 penguin) 



10 - ^ * 



6 8 10 12 14 16 18 

Running average of instantaneous luminosity (10^° cm'^ s'^) 

FIG. 3. Level 1, 2, and 3 trigger rates for the photon + 2 track trigger as a function of 
instantaneous luminosity in Run IB (open points). Total trigger rates for each stage are also 
shown (filled points). 
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pp^lr».^_d =/track impact parameter 

FIG. 4. The B decay vertex and relevant quantities on the plane transverse to the beam. For 
clarity, only the B momentum and one of its' charged daughters are shown. 
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6 8 10 
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-1 



6 8 10 

M(7K*K") (GeV/c') 



FIG. 5. Top: 7r+ invariant mass distribution for 5° — 7K*°(— K tt"^). There is one 



candidate. Bottom: 'jK'^K invariant mass distribution for ^(pi^ K^K ). There are no 

candidates seen. 
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FIG. 6. Invariant mass distributions of the K~'k^ combinations for B —>■ e~D^(—>- K~Tr~^)X, 
decays in the Run IB (top) and IC (bottom) data. The right-sign distributions (points) are for 
same charge electrons and kaons, as should be the case if they are both products of the real B 
decay chain, whereas in the wrong-sign distributions (histograms) the kaon has opposite charge 
to the electron. By fitting a Gaussian and a straight line to the right-sign distributions we find 
40.7±7.3 and 27.4±6.2 candidate B e~D°{-^ K-Tr+)X events in Runs IB and IC, respectively. 
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K (cm) 

FIG. 7. Photon conversion vertex density in the x — y plane in the 74 pb~^ of CDF Run IB 
inclusive electron data. The fine structure of the CDF tracking detectors can be clearly resolved. 
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FIG. 8. Photon conversion vertex density in ther — z plane in the 74 pb^^ of CDF Run IB 
inclusive electron data. The fine structure of the CDF tracking detectors can be clearly resolved. 
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M(eVr) = M(eV) + Mj/^ (GeV/c^) 

FIG. 9. Top: e'^e~K~iT'^ invariant mass distribution for — K*^{-^ K'''k'^)^{-^ e^e^) in 
the 74 pb~^ of CDF Run IB inclusive electron data. Bottom: corresponding e^e~K~ invariant 
mass distribution for the B~ — J/ijj{-^ e^e~)K~ reference decay. There are 28.0 it 5.8 events 
after background subtraction. 
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FIG. 10. Top: e''^e~K'^K~ invariant mass distribution for Bg — (^(— K'^K~)^{-^ e+e~) in 
the 74 pb"^ of CDF Run IB inclusive electron data. Bottom: corresponding e'^e''K~ invariant 
mass distribution for the B~ — > J/ip{-^ e^e~)K~ reference decay. There are 35.0 ± 6.4 events 
after background subtraction. 
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FIG. 11. Top: e'^e'~pTr'~ invariant mass distribution for A° — A(— f'7r~)7(— e+e~") in the 
74 pb~^ of CDF Run IB inclusive electron data. Bottom: corresponding e+e~i^~ invariant mass 
distribution for the B~ J /ip{-^ e^e~)K~ reference decay. There are 24.0 zt 5.3 events after 
background subtraction. 
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FIG. 12. Dielectron invariant mass distribution of the J/il) ^ e^e~ candidates in the 74 pb~^ 
of CDF Run IB inclusive electron data. The number of the J/ip ^ e^e~ events obtained by fitting 
the mass distribution to a function of 2 Gaussians and a polynomial is ~ 8000. 
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